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Abstract. The gluino- induced contribution to the decay b — > sj is investigated in 
supersymmetric frameworks with generic sources of flavour violation. It is emphasized 
that the operator basis of the standard model effective Hamiltonian is enlarged and 
that a suitable redefinition allows to realize the usual scheme of leading and next-to- 
leading logarithmic contributions when QCD corrections are included. The effects of 
the leading order QCD corrections on the inclusive branching ratio for b — ► sj are 
shown. Constraints on supersymmetric sources of flavour violation are derived. 



1. Introduction 

Like the usual weak interactions, also new physics interactions give rise to an 'effective' 
Hamiltonian operator at low energies (well below My/)- It only contains the 'active' 
particles, that is quarks, leptons, photons and gluons. The details of the heavy (new) 
particles and forces on them acting disappear and their properties manifest themselves 
in the form and the strength of the various operators in the effective Hamiltonian. 

In the standard model (SM), the only operator at tree level comes from H^-exchange: 

TUff ~ Grfa-yMffi-ffl) (1) 

where the q are the various quark flavours and factors like the CKM matrix elements 
have been dropped. The form and strength of the W couplings are embodied in Gp and 
the left-handedness of the quarks. Because of the various gluonic corrections, the number 
of active operators is much larger than just the above four-Fermi term. Dimensional 
arguments and renormalization group techniques allow to classify them systematically. 

If there are new forces, they will yield not only different couplings strengths, but 
also different forms for the interactions. This implies that there are in general more 
operators than the one in eq. (|1|) and also more than those obtained in the SM after 
QCD corrections have been included. 
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In supersymmetric models, there are many sources of flavour violations relevant to 
the decay b — > S7 . We shall consider here only those due to a flavour violating gluino- 
quark-squark vertex |J, whose impact on the decay b — > 37 was first analyzed in ref. |4[]. 
In general, the gluino contribution to this decay is the largest, unless it is assumed 
that the only source of flavour violation is encoded in the superpotential and that the 
deriving electroweak-scale flavour- violating parameters in the scalar sector are small. In 
this case, the gluino contribution can be neglected, unless tan (3 is large || . A complete 
calculation of the b — *> S7 rate, including corrections up to the next-to-leading order 
(NLO) )in QCD, exists for a specific class of models in which the gluino contribution 
can be neglected ||. Until recently, no study existed on the interplay between the 
gluino contribution to the decay b — > sj and QCD corrections, not even at the leading 
order (LO), in spite of the potentially important role that the flavour- violating vertices 
gluino-quark-squarks may have for model building (see for example Ref. 0). 

The gluino contribution to the decay b — > S7 exhibits two special features, when 
QCD corrections are implemented. First, gluinos couple to both left- and right-handed 
quarks. Thus, many more operators than in the SM, in which only left-handed quarks 
are coupled to the W boson, are expected. Secondly, the presence of the strong coupling 
a s in the gluino contribution and in the QCD corrections immediately raises the question 
of how to order the different powers of a s and whether to include them in the gluino- 
induced operators or in their Wilson coefficients. These two difficulties necessitate a 
dedicated study of the QCD corrections to this contribution. 

2. Ordering the QCD perturbation expansion and the effective Hamiltonian 

In the SM, rare B decays are induced by loops with W bosons; the coupling strength 
is always Gp. The corrections are due to gluons and other light particles and give 
rise to powers of the large logarithmic factor L = logim^jM^). In the decay b —>■ 57 
only loops with gluons contribute, and thus powers of L and a s are related as follows ||: 

LO: G F (a s L) N , (N = 0, 1, ...) 
NLO: G F a s (a s L) N . 

These terms are summed (over N) using renormalization group techniques. This results 
in an effective Hamiltonian 

4C 

nZf = , (2) 

where V t b and V ts are elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and 
H is the 'scale'. The Wilson coefficients Cj contain all dependence on the heavy degrees 
of freedom, whereas the operators Oi depend on light fields only. When all order QCD 
corrections are included, the Hamiltonian is independent of the scale \i. At a specific 
order in QCD, however, there is a residual dependence on \i. In practice, this scale must 
be chosen to be around m b . 
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The operators relevant to radiative B decays can be divided into two classes: four- 
fermion current-current operators and magnetic operators 

O m = ™&M (s^P R b) V, u , (3) 

where V^ u is the field strength of a photon or gluon. All these operators are of dimension 
six. 

It is by now well known that, in the SM, a consistent calculation for b — * s'j at LO 
(or NLO) precision requires three steps: 

1) a matching calculation of the full SM theory with the effective theory at the scale 
ji = nw to order oP s (or a l s ) for the Wilson coefficients, where \iw denotes a scale 
of order M w or m t ; 

2) a. renormalization group treatment of the Wilson coefficients using the anomalous- 
dimension matrix to order a), (or a 2 s )\ 

3) a, calculation of the operator matrix elements at the scale \i = ^ to order a° (or 
a]), where fib denotes a scale of order m b . 

Matters can be different in other B decays or when other contributions to b — > s^y 
are considered. An example is the decay b — > si I. The first large logarithm 
L = log(m^/M^) arises without the exchange of gluons. This possibility has no 
correspondence in the b — > 57 case. Consequently, in the case of b — ► s £ £, the decay 
amplitude is ordered according to GfL (a s L) N at the LO in QCD and Gfoi s L(oi s L) n 
at the NLO. To achieve technically the resummation of these terms, it is convenient to 
redefine certain operators and their Wilson coefficients as follows 0: 

ar = ^-ot, cr w = ^Q (< = 7,...,io). (4) 

This redefinition allows us to proceed according to the above three steps when calculating 
the amplitude of the decay b — > s££ @. In particular, the one-loop mixing of the 
operator O2 with the operator Og ew appears formally at 0(a s ). 

Including gluinos, we can now write the complete effective Hamiltonian as: 

H eff = HY ff + Hl ff , (5) 

where is the SM effective Hamiltonian in (0) and 7~l 9 e ff originates after integrating 
out squarks and gluinos. Note that 'mixed' diagrams, which contain, besides a W boson, 
also gluinos and squarks, give rise to a s corrections to the Wilson coefficients in H^ft 
(at the matching scale). Such contributions can be omitted in a LO calculation, but 
they have to be taken into account at the NLO level. As for Ti-lff, the aim is to resum 
the following terms: 

LO: a s (a s L) N , (AT = 0,1,...) 
NLO: a s a s (a s L) N . 

While H-eff i s unambiguous, it is a matter of convention whether the a s factors, peculiar 
of the gluino exchange, should be put into the definition of operators or into the Wilson 
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Figure 1. Diagram mediating the b — > 57 decay through gluino exchange and 
contributing to the operator Ojb,g- A contribution to the primed operator 0' 7b g is 
obtained by exchanging L <-* R. 



coefficients. It is convenient (and possible) to distribute the factors of a s between 
operators and Wilson coefficients in such a way that the first two of the three steps in 
the program for the SM calculation also apply to the gluino-induced contribution. This 
implies one factor of a] in the definition of the magnetic and chromomagnetic operators 
and a factor a 2 in the definition of the four-quark operators. With this convention, the 
matching calculation and the evolution down to the low scale /x& of the Wilson coefficients 
are organized exactly in the same way as in the SM. The anomalous-dimension matrix, 
indeed, has the canonical expansion in a s and starts with a term proportional to a\. 
The last of the three steps in the program of the SM calculation requires now an obvious 
modification: the calculation of the matrix elements has to be performed at order a s 
and a 2 at the LO and NLO precision. With this organization of QCD corrections, the 
SM Hamiltonian Ti^ff in eq. (||) and the gluino-induced one H^tt undergo separate 
renormalization, which facilitates all considerations. 

The list of such redefined operators is lengthy and given in ref . . For illustration 
only some representative ones are shown here: 

• magnetic operators, with chirality violation coming from the 6-quark mass: 

7b ,- g =eg 2 M fn b (fx) (saTnb) F„ v , (6) 

• magnetic operators in which the chirality- violating parameter is the gluino mass nig, 
included in the corresponding Wilson coefficients: 

O m =eg 2 s ^)(rsa^P R b)F llu , (7) 

• magnetic operators, with chirality violation signalled by the c-quark mass: 

7c , g = e g 2 s (fi) m c ( M ) (sa^P R b) F, v , (8) 

• four-quark operators with vector Lorentz structure: 

0\ W9 =g A M(s^P L b) (?fP L q) , (9) 

• four-quark operators with scalar and tensor Lorentz structure: 

O q 15 ,~ 9 =gt(v)(sP R b)(qP R q). (10) 

These last four-quark operators are induced by box diagrams only and through the 
exchange of two gluinos, whereas those with a vector structure are induced by box 
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and penguin diagrams. For each of these operators there are several similar ones, with 
different helicity or colour structure. Those with opposite helicity structure, obtained 
through the exchange L <-> R, are called hereafter "primed" operators. 

The four-quark operators in (|9]) and (|T0| ) are formally of higher order in the strong 
coupling than the magnetic and chromomagnetic operators. The scalar/tensor operators 
mix at one loop into the magnetic and chromomagnetic operators, giving rise to a large 
log L. Given this fact, the necessity of including 7c ~ g and 8c ,g in the operator basis 
becomes clear immediately. Due to these mixing effects, the scalar/tensor operators have 
to be included in a LO calculation for the decay amplitude. The remaining four-quark 
operators with vector structure 0\ lg -0\ ig (and the corresponding primed operators) do 
not mix at one loop neither into the magnetic and chromomagnetic operators nor into 
the scalar/tensor four-quark operators. Therefore, these vector four-quark operators 
become relevant only at the NLO precision. 



3. Wilson Coefficients at the Decay Scale 

At the electroweak scale, the Wilson coefficients of the various operators are obtained by 
perturbation theory; gluino exchanges such as in fig. 1 contribute. The renormalization 
program then allows to calculate the coefficients at the relevant scale m^. As already 
mentioned, the two terms Ti^t and T~L 9 e ff in the effective Hamiltonian @ undergo 
separate renormalization. The anomalous-dimension matrix of the SM operators 0\- 
(D$ and the evolution of the corresponding Wilson coefficients to the decay scale /if, are 
very well known and can be found in [|l(J. The evolution of the gluino- induced Wilson 
coefficients C it g from the matching scale nw down to the low-energy scale fi b is described 
by the renormalization group equation: 

= Q,sO) 7ji,s0) • (11) 
The usual perturbative expansion for the initial conditions of the Wilson coefficients, 

CtsiPw) 

Air 

as well as for the elements of r yj ig (fj J ), 

= iji,~ g + ina + , (13) 

is possible because of the choice of including appropriate powers of g s (fi) into the 
definition of the operators O^, as discussed previously. The anomalous- dimension 
matrix 7^^ is then a 112 x 112 matrix. However, the vectorlike character of QCD 
and dimensionality arguments reduce it to two identical 54 x 54 matrices. At LO, a 
further reduction to 8 x 8 matrices occurs. At the low scale, the LO expression for the 
Wilson coefficients of the dimension-six operators are 

C 7 b,g(Hb) = C 7Kg (n W ) + (77§ - 7?i) C %Kg (n W ) + R 7 b,g(Vb) , 

37 

C& b ~g{Vb) = V™ C 8brg (/j, w ) + R & b,g{Hb) ■ (14) 



CwM = <%fr w ) + Cl~(» w ) + , (12) 
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where the remainder functions R7b,g{f^b) and Rsb^i^b) contain the coefficients of the 
scalar/tensor operators and are small in this case. At the NLO, also the operators with 
vector Lorentz structure enter. 



4. Flavour violation and results 



While the main purpose of this talk was to stress the importance of including and 
systematically treating the new operators arising in beyond the SM physics, a few specific 
results should also be mentioned. Gluino exchange gives rise to flavour changes if the 
quark and squark mixing matrices are different (more precisely: do not commute). 
This happens when the so-called soft contributions to the squark mass matrix are 
substantially off-diagonal. Of course the values depend on the specific model; here 
we introduce phenomenologically the quantities (see ref. fill): 

r ( m d,£i)ij x ( m2 d,RR)ij C_L -\ /ie\ 

OLL,ij = 2 ' °RR,i3 = 2 ' K l T J) ( 15 J 

m~ m~ 

which are a measure for flavour violation if % ^ j. The various m? are the elements of 
the soft mass matrices. The parameters are varied, as usual, in their allowed ranges. 

Figs. 2 and 3 exemplify the results of this work. While the first shows how much 
systematic QCD corrections influence the branching ratio, fig. 3 illustrates its sensitivity 
to the soft matrix elements. 
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Figure 2. Gluino-induced branching ratio BR(B — > X s j) as a function of x = m|/m|, 
obtained when the only source of flavour violation is Slr,23 (see text), fixed to the value 
0.01, for Tfiq = 500 GcV. The solid line shows the branching ratio at the LO in QCD, 
for lib = 4.8 GeV and \iw = Mjy; the two dotted lines indicate the range of variation 
of the branching ratio when spans the interval 2.4-9.6 GeV. Also shown are the 
values of BR(_B — > X s j) when no QCD corrections are included and the explicit factor 
a s (/i) in the gluino-induced operators is evaluated at 4.8 GeV (dashed line) or at M\y 
(dot-dashed line). 




Figure 3. BR(B — > X s -f) vs. 5ll,23, when Sll,23 and Slr,33 are the only sources 
of chiral-flavour violation. The dependence on 5ll,23 is shown for different values of 
Slr.33- (solid line), 0.006 (short-dashed line), 0.01 (dot-dashed line) , 0.1 (long-dashed 
line). The value of x = m'~/m'~ is fixed to 0.3 and m q to 500 GeV. The vertical band 
indicates the experimental constraint. 



